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Free Motion of a Sphere in a Rotating Liquid parallel to the 

Axis of Rotation. 

By S. F. Grace, University of Liverpool. 
(Communicated by G. I. Taylor, F.K.S. Eeceived May 1, 1922.) 

§ 1. The motion is supposed to be a small disturbance from one of uniform 
rotation like a rigid body, small motion being defined to be such that the 
squares and products of velocity and vorticity components may be neglected 
in the expressions for acceleration. The system is supposed free from bodily 
forces, and is initially disturbed from relative rest by a motion suddenly 
communicated to the sphere. 

The pressure intensity of the liquid consists of two parts, one depending 
only on the distance from the axis, and the other on the disturbed motion. 
If at any instant the sphere is moving parallel to the axis of rotation, and 
the disturbed motion of the liquid is symmetrical with respect to a line 
through the centre of the sphere parallel to the axis of rotation, the motional 
part of the pressure intensity will also be symmetrical with respect to this 
line, and its resultant effect on the sphere will be to produce an acceleration 
parallel to the axis of rotation. Now the positional part of the pressure 
intensity would maintain any portion of the liquid in relative rest, and will 
consequently maintain the sphere in relative rest, providing the density of 
the sphere is equal to that of the liquid, which we shall suppose to be 
the case. 

The initial disturbed motion of the liquid will be irrotatiohal, since the 
effects of rotation take time to develop, and if the initial relative motion of 
the sphere be parallel to the axis of rotation, the initial disturbed motion of 
the liquid will have the symmetrical character above specified. 

It follows that the sphere will continue to move parallel to the axis of 
rotation, and the disturbed motion of the liquid will continue to be 
symmetrical with respect to the line of the sphere's relative motion, which 
in the sequel is called the axis of the sphere. 

The method of discussion will be to assume the particular type of motion 
indicated above, and then to show that all the conditions can be satisfied by 
making it a particular function of the time. This leads to a perfectly 
definite mathematical problem, the solution of which gives the following 
results : — 

The sphere oscillates about a point on its axis, the distance of which,, 
from the initial position of the centre of the sphere, is proportional to the 
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velocity of projection ; the amplitude of the oscillation tends rapidly to 
zero, but the period tends to a constant value which is half that of the 
rotation of the undisturbed liquid. 

With regard to the velocity of the liquid a general expression has been 
obtained for its value at any point, reducing to simple expressions in particular 
cases. These have been examined. 

On the axis of the sphere the velocity of the liquid is always less than 
the initial velocity of projection of the sphere and tends to zero at infinity. 
In the neighbourhood of the sphere the velocity oscillates about the zero value 
with a period which is ultimately constant and equal to that of the sphere's 
oscillation; the amplitude also tends to a value which is constant at a 
particular point, but there is an increasingly large gradient of velocity along 
the axis. 

On the equatorial plane of the sphere it appears that the disturbance is 
practically confined to the sphere's immediate neighbourhood. Away from 
the sphere the velocity of the liquid rapidly tends to zero, oscillating about 
the zero value with decreasing amplitude and a period which is ultimately 
constant at a particular point. At infinity this period becomes that of the 
ultimate oscillation of the sphere, but increases without limit as the sphere is 
approached. The velocity gradient over the plane also increases without 
limit. 

Finally the velocity of the liquid tangential to the sphere rapidly approaches 
a value which is constant in magnitude at any point, varying from zero at the 
poles to a maximum at the equator, the direction, however, changes with 
constant angular velocity which has a maximum value at the poles and is 
zero at the equator. Here again there is an unlimited increase in velocity 
gradient, in this case along a meridian of the sphere. 

However, since the gradients of the velocity of the liquid over the equatorial 
plane and of the transverse velocity along a meridian of the sphere ultimately 
increase without limit, parts of this solution are not applicable, for a stage is 
reached after which the components of vorticity cannot be considered as small. 
Consequently the physical assumptions underlying tlie mathematical treat- 
ment, viz., those of small motion, are violated ; in other words, the solution 
will only represent the true state of the liquid for a restricted time. 

The question of the ultimate physical state thus remains unanswered. 
J. Proudman* showed that a small steady disturbance was impossible, but 
then G-. I. Taylorf obtained a solution of the general equations of steady 

* " On the Motion of Solids in a Liquid Possessing Vorticity," ' Roy. Soc. Proc.,' A, 
vol. 92, 1916. 

t * Proc. Camb. Phil. Soc.,' vol. 20, Part HI. 
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symmetrical motion about a sphere, which does not satisfy the equations 
of small motion. Whether or not Taylor's solution represents the ultimate 
state in the problem of the present paper is a question for further 
investigation. 

The work was sugggested to me by Prof. J. Proudman, to whom I am 
greatly indebted for his interest and advice ; it is an extension of some of his 
own work (loc. ciL). It is proposed in a future paper to give the corresponding 
results when the sphere is projected at right angles to the axis of rotation, 
part of the initial work for which has already been effected. A combination 
of the two cases will give the small general motion of a sphere projected 
with given initial velocity. 

Mathematical Problem. 

§ 2. We shall use Cartesian axes Oxyz rotating about Oz with constant 
angular velocity o>, called the primary axes. Let %i> v, w denote the com- 
ponents of the velocity at x, y, z relative to these axes. The equation of 
continuity for a uniform liquid gives 

^ + ^ + ^ = o (1) 

dx dy dz 

and the boundary condition may be written 

lu + mv + nw = nY, (2) 

i, m, n being the direction cosines of the normal to the sphere and V the 
velocity of the sphere. These relations being kinematical are independent of 
the fact that the axes are rotating. 

The dynamical equations, in the case of small motion under no forces other 
than the pressure of the liquid, reduce to 

dt dx 

S + 2»*=-^ S (3 

dw = _3P 

dt tiz 



j 



where P =jp/p— ^co 2 (x 2j ty 2 ), (4) 

the letters p and p having their usual signification. 

Take a set of axes, having their origin always at the centre of the sphere 
(the secondary axes), and let these secondary axes be always parallel to the 
primary axes. 

In equations (l)-(4) the actual values of x, y y z are only required in the 
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term \ <o 2 {x 2 -\-y 2 ) in (4), since in 3/3#, d/dy, djdz we have t constant. Hence 
if x, y, z now refer to the secondary axes, the equations (1), (2), (3) remain, but 

P = p/ p — i a) 2 (x 2 + yo 2 l 

%o> y<$> % referring to the primary axes. 
The term ^ pco 2 (x 2 + y 2 ) in 

p = pT + 1 pco 2 (%o 2 + y 2 ) 

would maintain a liquid sphere in relative rest, i.e., would produce an accelera- 
tion a) 2 (xo 2 +yo 2 )i towards the axis of rotation. It will therefore produce 
an equal acceleration in the solid sphere, and we may take 

p = pV 

as the effective pressure intensity in producing the acceleration of the sphere 
other than centripetal. 

Now assume that the functions u, v, w, V, and P may be expressed as power 
series in time. This proves to be legitimate, since the results obtained show 
these series to be highly convergent for finite values of t. Let 
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P = 2o) X" P. - 2 ^ ); 

s = 



S i 



where -&*, ^, w a , c 5 , and P s are functions of x, y, z only, the initial state being, 
given by u , v , w , c . 

Substituting these expressions in the equations (1), (2), (3), we obtain 



dtt s dv s div s 
~r --s — r 



o, 



12 

V/J -I-J Ami J •»« ft 



3^ 3^ dz 

since the equation (1) is true for all time ; the boundary condition gives 

ht s -\-mv 8 -\-nio s = ?m? s , s = 0, 1, 2, ... , 

and from the dynamical equations 

u s+ i—v s = —cFs/dx^ 
Vs+i + Us = -o^s/dy 
w s+1 — —dP s /dz. 
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which remain valid for s = 1 if we replace dV-i/dz by — w Q . 
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The velocity of the liquid at any point may now be expressed in terms of 
the sequence P s , thus* 

u-~ tcojo (2g>£) — I'oYi (2ayt) 



co . co 



3/cte X F s y s+1 (2oot) — dfdy X P*7*+2(2a>0> 

s = s = 



v — %7o (2«0 + %7i (2©0 



00 CO 



= —d/dy £ 'P s y s +i(2a)t) + d/dx£ P,Y,+2(2©0, 

S = 5 = 

oo - 

w-wo = -3/a» X v s (2coty +i /(s+i) i 

S a 

CO 

where 7, (») = ^ ( — l) n s* +2n /(s + 2ri) I 

% = 

and we now concentrate our attention on the sequence P |f 
Inside the liquid we have, from (5) and (7), 

V *F S = -a»P,. a /&», s = 1, 2, 3 ... , 

and from (6) 

y 2 P = 0, 
while at the boundary 

3P ?)P 

5 = -rW^+i -%<?,-! — W -^=- 2 , S = 1, 2, 3, ... , 



&' 



^ ft 

cm account of the symmetrical character already referred to, with 

3P /3% = ~nci. 

The initial motion being irrotational, we may let P-i denote the velocity 
potential, so that 

u = — 9p_ 1 /aB, % = — 3P-i/3y, ^0 = — 9P-i/3«, 

and v 2 P-i = 0, 9P-i/ck = — w^. 

The value of c > the initial velocity of the sphere, determines the magni- 
tude of the subsequent motion; in the sequel this quantity, c , is denoted 
by W. 

The resultant pressure on the sphere is in the direction 0^ and equal to 

— /> JJ (P) r=a cos 0<#S 

taken over the sphere, where r, 0, <£ are the spherical polar co-ordinates of the 
element ^S referred to the centre of the sphere as origin, (9 = denoting the 
z axis and <j> = the plane xQz, a being the radius of the sphere. This 
becomes 

(2<oty 



20) P 2 e f 

5 = Si 



(P,) P=sa cos0rfS, 



* Proudnian, ?oc. tit. 
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and the equations of motion of the sphere reduce to 






bo that 
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c s+l = — JJ(P,) r = o cos0dS/($7ra 3 ), s = 0, 1, 2, ... , (8) 

iving the functions c s in terms of the sequence P*. 
Assuming the liquid to be at relative rest at infinity, the boundary 
conditions may now be written. 
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The function P has to satisfy these two relations and vanish at infinity ; 
this is impossible unless P is zero. Generally 
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^2s = 0. 

Hence we require a sequence of functions P&-1 (s = 1, 2, ....), which 
vanish at infinity and satisfy 

v 2 p 2s -i - -arw^ 

inside the liquid and 
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The Solution. 

§ 3. After an investigation which would take too much space to reproduce, 
the required sequence was found to be as follows : — 

P 2s -i = Wa { (ajr) 2 Q 2 , 2s (cos 0) + (a/r) 4 Q 4 , 2s (cos 0) + ... 

. . . + (a/r) 2 ™ Q 2m 2s (cos 0)+... + (a/r) 2s+2 Q 2s+2 , 2s (cos 0) }, (13) 
where, putting cos = fi, 

q„ „ oo = ( - ir ■ 1 • ^i) ; - ( ;i 2 ; - + ^ , * a - ^- x 



_2m + l L 3 " ' 3 

,/ iv »-i2s+l 2s + 3 28+2ot-3 8m - a 



5 2m -1 



2s+l,„ r , s , 2s + 1 2s+3 

3 " - 3 



l_f£±i»C 1/A 2 + ^ii± . fi±i»c 8 M*' 



,/ iy»2s + l 2s + 3 2s+2m-l 2m \ 1 n „ 
••• + (-1) -3-.— g-... 2m + 1 /* J-j, (14) 

for s> and m = 1, 2, ... (s+l). A more compact form is given by 

Q^,.^) = (-1)— | . —^J—j— x 

(^r 1 !^^-^- 1 -^ 1 ^-^}' ( i5 > 

and it will be noticed that 

Q*+2,2.(/*) = (-1)^P2.+ i(/a)/(2s + 3), 

where, in this case, P 2s +i (/a) denotes Legendre's function. 

We shall now simply verify that these formulae satisfy all the conditions. 

Verification. 
| 4. Considering first the relation (9), the expression (13) satisfies this 
provided 

V 2 {(a/r) 2 ™Q 2m , 2 s} = -d 2 /dz 2 {(a/r) 2 ™Q 2m , 2s - 2 }, s = 2, 3, ... , 

^.(l^^ ) 3^^_ 2 ^^|^4-2m(2m~ljQ 2m , 2s 

+ 2m{2(m + l)/i 2 — 1} Q2m,2 5 -2 

This relation is satisfied when m = s+1 since Q 2s+2 , 2 * is Legendre's 
function of order (2s + 1), apart from a factor, and Q^+^aj-s is zero. 

Otherwise let 

n — ( ly-Bt^ 1 .3 . 5 ... (2s— 1) 

Qa», a. - ( - 1) § 2*(m-l)!(s-m+l)! ?2m ' 2s ' 
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then we have to show that 
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and 



Hence 

( 1 -V) ^|^- J - 2ya S -^ s + 2m (2m - 1) q 2m , 2s 

= (1 - M 2) S -»»-i [ {2s (2s-4m + 1) fi 2 + 4m 2 - 2s} * 2s 
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and the coefficient of /^ +1 In (17) is 
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and it is easily verified that the coefficients of p, /J, /^ 2m+1 , yu 2wl+3 in (16) and 



17) are also identical 
Examining now the case where s = 1, we have 
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v ; ^ 2 ^ v H<) ( 2 s+l)(2s + 3) (2s +1) (2s + 3). 

= (-l) , W[i.^(l- / *»)-3/{(2* + l)(2« + 3)}]. 

H 2 



100 Mr. S. F. Grace. Free Motion of a Sphere in a 

Substituting now in equation (10), we have 
(-iyW[|^-i(l- / . 2 )-3 / i/{(2 S + l)(2 S +3)}] 

3W . , „,._„ 3W 



H(-«) 



s—1 



+ (-1) 



s-2 



(2s-l)(2s+l)J 



(2s +1) (2s + 3) 

-M-rw {|/* a - s (i-M s ) 3 



for ^ = 2, 3, ... and 

aPi/8r> =0 = -W{f M (l-^-*/*> 

cos# 



(2s-l)(2s+l).. 



W{j^-f^} 



|-7ra 3 



(Pi) r = « cos 6 dS — W cos — cos # 



| ^JWa(-Ycos0 



.3# 1/ W j _ir=^ 

It follows that the sequence P2«-i, given in (13), satisfies the conditions of 
the problem, and therefore determines the motion, 

Convergence. 

§ 5. Before proceeding to the results obtained from this sequence, we shall 
consider the convergence of the series for u, v, iv, and P, in terms of P^-i. 
From the expression for Q 2m > 2s (/a), given in (14), it is seen that, for m^s 

Q a »,*00l<2 • 2* (m- 1)! (5-^+1)1 X 

[2(2m + 3)(2m + 5) ... (2s + 2m- 3) {l + m " 1 Ci + 8| - 1 C 3 + ...} 

+ (2m + 3)(2m + 5) ... (2s + 2m-l) {l + m Ci + m G 2 + ...}] 
<-§ (2m + 3)(2m + 5) ... (2g + 2m-3) „ g . 

3 o^.j (s + m)(s + m~ l)(g + m— 2) ... (»i-f 2) 
"^2 (m— l)!(s-m+l)! 



<o 



€> 2 m ~ 1 



2(m-l)! 



2s(2$-l)(2s-2) ... (s + 2). 



Hence 



Qa»i, 2s (/*) 



<C^ o . 



/i 



5 



(25-1)1 "^ * (m-1)! (s-fl)l 
and | Qss+2, 2 S 0&) | =s 3/ (4s 4- 6). 

Consequently 



m < s 



1 



•* 25-1 i ^ j^ yf 



9m~l 



< 



24 



a#d 



W« (2s— 1) I ^ s ! m== i (m— 1) I 'si 

(X) 

(2©0 a> " 1 /(2s-l)! 



P| - 2a> JjT jX"2s-) 

5 = 1 



oo 

< $ 48Wao> (2aty- 1 /s ! 



S = 1 



Rotating Liquid parallel to the Axis of Rotation. 101 

i.e. the series for P in terms of the sequence P25-1 is absolutely convergent for 
finite values of t. 

Similarly it may be shown that 

I Gradient of P25-1 1 < (4s + 3) | P25-1 1 , 

and that the series for u, v, w y in terms of P2 S -i, are absolutely convergent for 
finite values of t. 

The Motion of the Sphere. 

§ 6. The velocity of the centre of the sphere is given by 



OD 



V = X c, (2<oty/s ! 

and from (19) 

c 2s = (-l)»3W/{(2* + l)(2s + 3)h 

c 2s +i = J 

Hence 

•tr ^r* / -t\$ o W \2cot ) s 

^ s = o^ ' (2s+l)(2s+3) • l2i)T 

— SW J" sm 2ai cos 2a)< "\ 
~~ L (2<»0 3 ~ (2«0 a J * 

This may be expressed in a neater form by using the notation 

S B (a>) = a;» +1 (-1 . iLY 1 s i!Lf 

\ <JG (mO 1 Jj 

so that V = 3W . Si (2o)t)/(2cotf. 

The values of V/W for 2o>£ = 0, 1, 2 } ... 14 have been found, and from them 
the Diagram I constructed, which shows the ratio V/W plotted against that 




o 



0-5 

Values of ^/f 

Diagram I. 
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of t/T, where T = 27r/« is the time of a revolution of the undisturbed 
liquid. 

It is to be observed that the value of V is oscillatory, but tends rapidly to 
zero ; the sphere oscillates about its final position with continually decreasing 
amplitude, the period of oscillation tending to a constant value equal to ^T ; 
after one revolution of the liquid the sphere is virtually at rest. 

The distance between the initial and final positions of the sphere is given by 



'CO 



roo 



Vdt = 







3W 



o 



{ 



sin 2cot cos 2a)t 
{2wtf ~~ ' (2cot) 2 



at ~ 2 ' 2 * 2o> ""16 



The quantity WT is the distance the sphere would travel during the time T 
with its initial velocity of projection, consequently the sphere comes virtually 
to rest after a time T, at a point 3/16 of this distance from the point of 
projection. 

Velocity of the Liquid along the Axis of the Sphere. 



§ 7. Along the axis of the sphere we have fi 
case /a = 1, we derive from (15) 



+ 1, and taking first the 



Q 



2nu 2s 



0, 



m < s, 



with 



so that 



Q 2s , 2s = (-l)'- I 3/(4«+2), Q 2s+2;2s = (_l)'3/(4« + 6), 



U-x-C-l)— ^-j 2JT3 



'a\ 2s 1 



\T J JiS "T* -*- 



and on the axis of the sphere 



dz Br 



v ; 2 XW 2^ + 3 \rj 2^4-1 



From equations (6) we have 



w s = — d~P s ~i/dz, s = 0, 1, 2, ... , 



hence it follows that 



Since 







8P 



2s- 1 



CZ 



= (-!)• 



3W f /a\ 2s ^2s + 2 M 2 * +1 2s 



2 L 



a\ 

x r I 



2s -f 3 \r) 2s+l 



CO 



w 



Jjf w s {2cot) s /sl 



5>0. 
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we have, on substituting the expressions obtained for the coefficients w % 



QQ 



w 



s = 

~~ 2 \r 



2 1 W 2s + 3 W 2s + 1 J (2s) ! 

f_ 1 v 2s 4- 2 (2a(otfr) 2s 
( j 2s + 3 (2s) I 



s = 



2s (2a<ot/r) 2s 



-t--« --S^/ ^ 2s + 1 (2s)! 



2 r 



3W /«\* 



COS 



<9a,*Hr\ sin (2a<x>tfr) o cos (2aa)t/r) 



, sin (2ao)tfr) \ SW a f /0 . # N sin (2acot / r)' 

4. 2 — ^ — : ; ( o < ^ — -w - -< cos (zaootfr)— *— ^ 



(2acot/rf 



2 r L 



(2acot/r) 



3W fa 



9 



r 



«r+2* 



r (2a)tf) 



1 1 q I 2awt \ 

vtfj \ r J 



(Ajl \ 



which, on r = a, agrees with the expression obtained for the velocity of the 
sphere. 

The case /& = — - 1 gives the same expression for w, and, on the axis of the 

sphere, 

to = v = 0. 

Perhaps the clearest way to discuss the function w is to examine the 
variation with r for constant values of t A table has been constructed and 
from it the Diagram II formed, in which the ratio w/W is plotted against 
that of rja ; the curves are numbered 0, 1, 2, . ... 10, and represent tofW at 




Values of Va 
Diagram II. 



104 Mr. S. F. Grace. Free Motion of a Sphere in a 

points on the axis of the sphere for 2mt = 0, 1, 2, . ... 10 respectively. It 
may be noted that 2cot = 10 corresponds roughly to t = 0*8T. 

A further Diagram III has been drawn showing the velocities at times 
corresponding to T, 2T, 3T, respectively. 

It may be noted that the function Si (2aa>t/r) oscillates between bounds in 




Iff 



J*^lwl.M^<rf>^ | IS B '*>.»*»^i I hi. 




z> 



4 



5 



Values of T / a 
Diagram III. 



the neighbourhood of ±1. Eeferring now to the Diagrams II and III and 
to the expression for w given in (21), it is seen that w is always less than W. 
As r increases, w is ultimately positive and tends to zero. At intermediate 
points the velocity may be positive or negative, the disturbance being a 
maximum in the neighbourhood of r = Via. Here it is to be noted, however, 
that as t increases the velocity gradient along the axis becomes large, due to 
the fact that in dw/dr there is a factor 2wt. This point is clearly illustrated 
in Diagram III, which indicates the state of the liquid as t becomes large ; 
the changes of velocity are observed to be contained within upper and lower 
bounds, which are ±|W{ff/r-(ft/f) 3 }, since for t large we have approximately 



W 



2 \r W v \ r }' 



Velocity of the Liquid on the Equatorial Plane of the Sphere, 
§ 8. The equatorial plane of the sphere is the plane defined by fi 
Here we have 

: V = 0, 



U 



so that the velocity is wholly parallel to the axis of rotation. 

\ dz /> = T \ d/JL /fi = 



We have 



0. 
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and from (14) 

vQ2m,2s\ -/_iy-ro^ 1.3.5. ... (2s — 1) f 2^ 

'' - -^ } 2 # 2«(m-l)"- - ' 1Xf ^ 7r ~' 

- m+1 3 2m— 1 1 . 3 . 5 . ... (2s- 1) 



^ /m = o ' ~ 7 2 " 2 s (m — 1)! (s — m + 1)! l_2m-fl 

= (-iy 



IS- 5, m > 



so that 



2 * 2m +1 * 2*(m-l)!(s-m + l)! 



,_ 1V 3W 1.3.5....(2s -l) 

&7m = o"" 1 ; ~ 2^ X 

1 1_ M 3 _3 1 /g\» 5 1 (a* 
\ 3 # s I \r) 5 (s-1) ! 1 ! \r) 7 (s-2) ! 2 ! ^ 



..-(-!)*» 



1 



2m— 1 _^ 

2m + 1 (s— -m+l)l (m— -1) ! W 



' a \2m+l 



— 






— Z' 1 \s * * * * • • v ^ """"" -^ / 

-i- ) -g- • 27T71 



x 



-^/-•^(^•^(f 



<_1)..CV, ?2^4^r , + ... + (-l). 2£±ifer r 



2m + 1 W 



25 + 3 W 



The quantity within the brackets, on putting a/r = X, becomes 

1 tf-'d | \ 5 -f *C 2 5 \7_ _ + («iy |£+1 x2f+8 



(22) 



3 



5 



7 



25 + 3 



Hence 

\ 9^ /,* = 



\ 2 (l~\ 2 ) s d\-2s [ X^l-X 2 )^ 1 ^ 

o Jo 

X, s (l-A, 2 )*-2 f \ a (l-\»)»<Z\. 

Jo 



23) 



, 1V 3W 1.3.5. ... (25-1) 

- (~i) -£- ¥T7i 



x 



\3(l-\»y-2 



'A. 






so that, on the equatorial plane, 

CO 

w = % w s (2cot) s /s I 



\ 2 (i-\*yd\ k 



s = 



j WX 8 -^" ( _ l). 3W . 1-3-5 (2g-l) 

^ A , S I 



X 



s = l 



{x»(l-X-).-2^(l-X-).dx}^. 
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Further, 

V* ( i v - : 2. ' ^ •" (^ s ~"^) (\ w (2(ot) * 
** l ; 2*.s! l M (2s) ! 



,5= 1 



Jo {2&>£ ^/(1 — X 2 ) } being the Bessel function of zero order. Hence we have 



5 = 1 



w - _AW\»-JW 



\ 3 Jo{2^. v /(l-A 2 )} 

-X 3 -2 f A [X 2 J {2 < «if v / ( 1 - ;v - 2 )}-^]^ 



_Ji W 
2 ,v 



X 3 Jo {2a>* v/(l-X») }- 2 x2J o {2ft>V( 1 -^ 2 ) } <**• 



(24) 



This is the most compact expression which has been obtained for the 
velocity at points on the equatorial plane ; it' is convenient in that it enables 
us to obtain an approximation to the value of w. 

When X is small, the value of the integral is approximately -J-X 3 J (2cot),. 
provided 2(ot,j$ not large, which gives w in the form 

~i"WX s J (2«0. 
To approximate to the value of w for i large, we use the well-known fact 
that, for x large, J (x) is approximately 



2 \ / 7T 

— cos [so — — 
Sirxj \ 4 



so that, provided ^/(l - X 2 ) is finite and 2a>t sufficiently large, w becomes 
3WA 3 / f 2 1 f„ . ,„ .,i irl 




l^ 5 v /a->.') ; C08 l 2 " V(1 " x ' ) "iJ 

+ 3W |> V / {-^j-^} eos {2a,V(l-X 2 )-| I- «*, 



r v , 

JO V L7T 

and, on examination, we see that ultimately the integral is of a higher order 
of small quantities, so that w is approximately 
3WX 3 



2 




Wa^W 2 "'^ 1 -^-!^ 



It follows, therefore, that there is ultimately an oscillation about the 
equatorial plane, the amplitude at any point continually decreasing and 
tending to zero, the period, however, tending to the constant value 

2fi>v/(l-X T ) ~ 2 ' ^/{l-(a/ryy 

Again, since dw/dr contains a factor (2(ot)i t the velocity gradient ulti- 
mately increases without limit. 
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To examine further the equatorial velocity, we put the expression (24) 
into a power series in X. To obtain this form we return to the expression 
for (3P2s-i/3^) M = o given by (22) and consider the coefficient of \ 2m+3 for 
m > in 

'aPa.-i\ (2cot) 2s 



S* 



fi = 1 



dz / m = o 2s ! 



It is 



( — l) m 



2m + 1 3 W 



CO 



2m 4-3 
2m + 1 



2 



2 (-iy*c, 



(2o»t) 



2s 



m 2 2s {s !)a 

00 

2« + 8-2»«n--2-< 2flrf >"^ (-!>• 



s = m 

1 3W 



(2o)/) 2 " + 



W 



3 W 2m + 1 



2 2m + 3 2 m m ! 
and the coefficient of X 3 is 



(20)0^ Jm (2c0t), 



-JW{1-J (2»0}. 

Hence w is given by 

|W |i\ 3 Jo(2^)+| . ix 6 2fi>*J 1 (2a>0 + | • 2^4 



X 7 (2a)0 3 J2(2to^)+... 



2m + 1 



X*»+3(2©0 w J»<2a>*)+... 



2m -f 3 2.4.6 ... 2m 

From this expression we are able to construct a Table giving the values 
of w for 2cdt = 0, 1, 2, ... 10, and to graph the function w for constant 
values of t (Diagram IV). The Diagrams II, III, IV are all drawn to the 




Values of T / a 
Diagram IV. 
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same scale for purposes of comparison, but, in the last case, only the curves 
corresponding to 2cot = 0, 2, 4, 6, 8 are shown, since the variation is so small 

An examination of Diagram IV confirms the results already given ; it also 
indicates that the disturbance on the equatorial plane is practically confined 
to the liquid in the immediate neighbourhood of the sphere. 

On the sphere itself, putting \ = 1 in (23), 

9P*-i \ = / iyh-1 3W 1 . 3 . 5 . ... (2s- 1) 2 1 2s 2s-2 2 

H5/m = o K } 2 2 s . si ' 2s + 3 " 2s + l ' 2*-l'"3 

= (-l) 5+1 3W/{(25 + l)(25 + 3)}, 
so that, on the equator of the sphere, 

w- 2 dvv -4^^ -U (2« + l)(& + 3) ; (2s)! 
- _$ W + 3W . Si (2cot)/(2coty } 
which rapidly tends to the limiting value — f W. 

Velocity of the Liquid at Points on the Sphere. 

§ 9. To discuss the velocity of the liquid at points on the surface of the 
sphere, we transform the dynamical equations (3) into spherical polar 
co-ordinates, and obtain 

— —2wv cos# = — ^, 
at rdv 

§£ + 2ww' sin 8 + 2cm' cos = 0, (25) 

vt 

dw' . 3P 

-* 2<ov' sin = -—7T-, (26) 

3£ 3r 

where w' coincides with w for = and — %' coincides with w on the 
equator. 

On the sphere, we have 

z / X ^r/3 ow a rsin2<w£ cos2<y£ 

<« V-. = V cos = 3W cos ^—-j.-^-^ 

and from (20) 

3P2s~l\ __ / 1 \. W f3 _2.,-.] /l • 2 a 3 COS 



'r~ a 



(-lywllco^-^rin^- ^ } *>0, 



t 
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so that 



oo 

2ct> J5f 



'3ft 



(2wty 



=zO\dr Jr~a Si 



00 



- Za,2, s = 1 { i)W^ 2 cos (fan tf (2s +1) (2s + 3) J (2s-l)! 



= -2«W 






~ sin 2 # sin (2g>£ cos 6) 

4- 3 cos # 



{ 



sin 2^ 3 cos 2o)£ _ 3 sin 2g)£ "\ 
(2mtf (2wtf (2wtf J „ * 



Hence from (26) 

('y'^^a = — -f W sin sin (2a>£ cos #), 
and from (25) 

(~Or=a = -f Wsin<9eos(2^cos60 + 3Wsin<9 . Si (2a>t)/(2a>t) 2 . (27) 

The expression for (w')r = a has already been treated; it is simply the 
radial component of the velocity of the sphere. 




Values of 

Diagram V. 
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Considering the expression for (v') r = a, a Table may be formed giving its 
value for between and \ir at times corresponding to 2cot = 0, 1, 2, ... 10 
respectively, and from it the Diagram V constructed. The values for from 
\ 7r to 7r are those from to \ it, repeated in the reverse order with change 
of sign. 

It is observed that the transverse velocity, (v') r = a> oscillates between 
bounds + § W sin 0, which are independent of the time. 

The velocity along a meridian, coinciding with w on the equator, is given 
by (27). The variation of this quantity is indicated by the Diagram VI, 
which is constructed similarly to V. The values for from Jtt to it are 
those from to -Jtt in the reverse order. 

The value of (— w') r = a oscillates between bounds in the neighbourhood 
of ±f-Wsin-0, approaching these bounds more nearly as t increases; the 




Diagram VI. 
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irregularity is due to the term 3Wsin# . Sj (2a)t)/(2a>t) 2 } which rapidly tends 
to zero. 

Hence, for t large, the velocity of the liquid tangential to the sphere at a 
point of co-latitude #, has components 

— f W sin cos (2cot cos 0), — f W sin sin (2oot cos 0), 

which combine into -|Wsin in a direction making an angle a, measured in 
the negative direction from the meridian (0 increasing) through the point 
where a = 2cot cos 0. 

Hence, ultimately, the velocity of the liquid tangential to the sphere at 
any point is constant in magnitude and equal to f W sin 0, but changes its 
direction with constant angular speed 2m cos } which is a maximum at the 
poles and zero at the equator. 

A point to notice, however, is the fact that (dv' /d0) r= = a contains the 
factor 2cot, and therefore ultimately increases without limit. This means 
that we reach a stage beyond which the components of vorticity can no 
longer be regarded as small. A similar deduction is obtained from dvj/dr 
over the equatorial plane, which contains the factor (2a)ff. A physical 
restriction to the solution of the problem is thus introduced. 



